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Abstract . The multivariable Instrumental Variable/Approximate Maximum 
Likelihood (IV/AML) method of recursive time-series analysis is used to 
identify the multivariable (four inputs-three outputs) dynamics of the 
Pratt and Whitney FlOO engine. A detailed non-linear engine simulation is 
used to determine linear enoine model structures and parameters at an 
operating point using open loop data. Also, the IV/AML method is used in 
a direct identification mode to identify models from actual closed loop 
enoine test data. Models identified from simulated and test data are 
compared to determine a final model structure and parameterization that 
can predict engine response for a wide class of inputs. The ability of 
the IV/AML algorithm to identify useful dynamic models from engine test 
data is assessed. 

Keywvirds. Multivariable identification; Jet engine; Closed loop data; 
TnTTnjmental variable/approximatn maximum likelihood 


INTRODUCTION 

A tvoical engine contral design cycle 
consists of developing a dynamic engine 
simulation f?-om steadv-state component per- 
formance data, designing a control based 
upon this simulation, and then testing and 
modifying the control in an engine test 
cell to meet performance reguirements. 

This design cycle has been successful for 
state-of-the-art engines. However, for 
more advanced multivariable engines that 
exhibit strong variable interactions, this 
procedure will result in substantial trial 
and error modification of the conLrol 
during the testing phase. One method to 
automate the design process and reduce con- 
trol modification testing and development 
cost would be to identify accurate dynamic 
models directly from the closed loop test 
data. These identified models would then 
he used in conjunction with a synthesis 
procedure to systematically refine the con- 
trol. Recent advances in closed loop 
idertif iabi 1 itv (Ref. 11 present a methcd- 
olooy for this direct identification of 
enoine model dynamics from closed loop test 
data. This paper describes the application 
of the IV/AML identification method 
(Ref. 21 to simulated and actual closed 
loop FIDO engine data (Ref. 3). This study 
was undertaken to determine if useful 
dynamic engine models could be identified 
directly 'rom closed loop engine test data. 

Recently, some gas turbine models have been 
identified. In Refs. A and 5 a single 
input state-space model was obtained by a 


time series method from simulated input- 
output data. In Ref. 6 a nonlinear filter- 
ing technigue was used to estimate system 
parameters from multiple input simulated 
engine data. In Ref. 8 a two input engine 
model was determined from open loop data 
analysis using the "Method of Models." 

Later the model form of Ref. 8 was applied 
to condition monitoring for a single input 
enoine (Ref. 9). In Ref. 10 a single input 
enoine model was determined from closed 
loop flight data using a maximum likelihood 
parameter search of dynamic engine simula- 
tion parameters. It was assumed that there 
was no process noise. In Ref. 11 a two 
input automotive oas turbine model was 
identified from open loop enoine freguency 
and step response data. The two data types 
were used in concert to fix the model 
siiucture a’d to fit the model parameters. 

In this report the IV/AML method is applied 
to both simulated and actual closed loop 
test multiple input data. The IV/AML 
method is an output error identification 
method and was implemented in a combined 
iterative/recursive form. As vet, the 
IV/AML meth(id has not seen wide application 
to many physical processes, and, until now, 
the method has not been applied either to 
oas turbine engine data or to closed loop 
data. The test data studied in this report 
contains both measurement and process 
noise. The available closed loop engine 
test oat? records are each comprised of 
onlv 200 sample points. Since this is a 
relatively low number of sample points per 
operating record, the IV/AML method was 



ORIGINAL PAGE IS 
OF POOR QUAUTY 


2 


^plprtpd for ii';p hprausp Monte Carlo tests 
have <;hown the method to exhibit reasonable 
converqence for a small number of samples 
(Rpf. 2). 

This report beqins with a brief description 
of the IV/AML method and its application to 
enqine model identification. Next, an 
enqine model structure is developed and 
applied to simulation and test data. 

THE IV/AML IDENTIFICATION 
METHOD 

The multivariable IV/AML method is the time 
series analysis tool used in this study. 

Monte Carlo simulation has shown the method 
to be an asymptot ical 1 v efficient identifi- 
cation tool. A complete descripton of the 
Refined IV/AML method 1 -. qiyen in Ref. 2 . 
Briefly, the model is written as: 

Al;'^ lX|^ « )U|^ (la) 

C( 7 '’)C^ e D( 7 '^)p^ (lb) 

where i ^ is the backward shift operator. 

"v ’’u 

The vectors *^1^**^ known 

ard represent the noisy output and deter- 
ministic input, respectively. The vector X|^ 
is tt'e deterministic part of ,V|^ and 
modris the effects of unmeasureable disturb- 
,’nces ant' measurement noise. The vector 
n 

e, fR '' sat isf ies: 
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The polynomial matrices A. B, C and D 
are det ined as: 
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where A^, B^, and Oi are real 
matrices. The algorithmic structure of the 
iterative/recursive version of the IV/AML 
method is given in Fig. 1 . Here the auxil- 
iary system model is updated iteratively 
after a complete passage through the data 
during which the model is maintained at the 
previous iteration value, 8. . The noise * 

°T 

model, however, is updated recursively. 

Model updates are determined by the IV 
algorithm and the AML alqorithm blocks of 
Fig. 1 . These two algorithms have the same 
basic structure and represent iterative 
solutions to the necessary conditions for 
maximizing a log likelihood function of the 
observations. The algorithms basically 
update the model parameters as a function 
of output error and estimated parameter 
uncertainty. Prefiltering, which is used 
in a "refined" IV/AML procedure, can yield 
estimates with improved statistical effi- 
ciency. Prefilterinq was not used in this 
study, however, due to the additional com- 
putational complexity and the satisfactory 
operation of the "standard" IV/AML approach. 
The noise covariance was estimated by: 

5k ■ Vi • r ['k'l - 3k-i] 

When applying the IV/AML method of Fig. 1 
to the enqine data, usually six iterative 
passes through the data were taken. 

Initial parameter values for the system and 
for the noise covariance were determined 
from open loop simulation data. Values for 
the initial asymptotic parameter covariance 
matrices, * for the parameters of 
and R for the e^, were selected as: 

t„ . .1 

( 5 ) 


where u ■ up in general and u is 
selected as approximately 100 to 1000 times 
larger than the diagonal elements of Q, 
Once estimates of #3, e^, and Q were 
obtained in one test, they were often used 
as initial values in subsequent tests. 
However, during each test, Q was normally 
held constant during the iterations 1-4 and 
6, and updated during iteration 5 . Eigen- 
values of A, P, and D were calculated at 
the start of each iteration to nsure 
stable alqorithm operation. Also, compari- 
son of eigenvalues from one iteration to 
the next is a good indicator of algorithm 
performance. For example, the irovement of 
one eigenvalue of % toward the stability 
bound indicates that the alqorithm is 
"getting lost" or that an incorrect model 
structure has been chosen. Also, norms 
P|j and R|j were calculated at each 
iteration as a measure of converqence. 
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ENGINE MODEL 

The Pr'tt and Whitney FlOO engine (Ref. 3) 
is a twin-spool low-bypass ratio after- 
burning turbofan. Four controlled vari- 
ables are considered: main fuel flow (WF) , 

exhaust nozzle area (AJl , compressor (fan) 
inlet variable guide vanes (CIVV) , and the 
rear compressor variable guide vanes 
(RCVV) . Three output variables are con- 
sidered: engine fan speed (Nl) , engine 

compressor speed (N2) , and augmentor 
entrance pressure (PT6) . The engine speeds 
are indicative of the dynamic response of 
the engine while PT6 is closely related to 
engine thrust. Globally the engine is 
modeled as 


X • f (x.u.ALT.MN) 

( 6 ) 

y s g(x,u,ALT,MN) 

where x is the state vector, u is the 
control vector and y is the output vec- 
tor. Engine operation is also dependent 
upon environmental variables altitude (ALT) 
and Mach number (MN) . An engine operating 
point is defined as 


flXss^^ss.ALT.MN) - 0 

q(x .u .ALT.MN) « y 
ss ss •'ss 


(7) 


A thirc order behavioral model relating the 
engine outputs to the primary control vari- 
ables WF ana AJ was developed in Ref. 12 
and is given as 
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This model represents linearized behavior 
in a small rqion about an operating point. 
Including CIVV and RCVV and writing in the _ 
form of (11 . the behavioral model becomes 

♦ A^z'^^x^ . BjZ'^^ (9) 

where x • (Nl. N2. PT6) and u ■ (WF, AJ, 
CIVV. RCVV). 

Si mu lation Application 

The IV/AML method using the model of (9) 
was initially applied to open loop SISC 
digitally simulated engine oata. These 
data were obtained for an ALT > 10 000 ft, 
MN « 0.9 intermediate power operating 
point. Two data sets were generated in the 
following manner. All controls were held 


constant at their steady-state operating 
point values except for WF in the first 
case: 


WF - WFjjd ♦ 0.03 Y^) (10) 

and AJ in the second case 

AJ - AJjjd ♦ 0.03 v^) (11) 


where y^^ is Gaussian noise with 

=K} ■ » 
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(12) 


Equations (10) and (11) represent *3 per- 
cent variations about the nominal operating 
point to Insure linear operation. In both 
cases the sampling rate was T - 0.01 sec. 
From these two data sets an initial esti- 
mate of Aj was obtained as 



0 

-0.96 

0 



(13) 


This corresponds to continuous eigenvalues 
of -16, -4, and -90. This information was 
used as a starting point in MIMO identifi- 
cation test. 


A random input emphasizes identification at 
all frequencies. However, the frequency 
range of engine operation is practically 
limited by actuator and engine hardware to 
about 20 radians/sec. Thus a set of inputs 
was chosen, somewhat arbitrarily, for a 
MiMO test that empahsized this more realis- 
tic frequency range. Here 


U|^ - Ujjd ♦ 0.03 aU|^) 
id 

‘(sin(O.lkT) ♦ cos(5kT ♦ 2) 
(sin(0.2kT) ♦ cos(7kT ♦ 4) 

au. > 

“ (sin(0.3kT) ♦ cos(5kT) 

(cos(0.4kT - 4) - cos(lOkT) 


(14) 


(15) 


Also the sampling interval was selected to 
be T ■ 0.05 sec to correspono to sub- 
sequent tests which used actual test data 
recorded at T ■ 0.05 sec. Four hundred 
data points (K ■ 4001 for each test were 
recorded. The IV/AML method was applied 
and the results are given as model 1 in 
Table 1. The average percent error per 
point was calculated from the outputs as 
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Error. > 


(16) 


and is given in Table 1 also. For model 1 
each output has less than 0.6 percent error 
per point on the average which is a good 
fit of the data. An additional element in 
Aj was found to be required to satis- 
factorily model PT6. 


The noise model (see Fig. 1) was found to 
be close to the auxiliary model. Thus, the 
number of parameters could be reduced by 
constraining Cj = Aj . However, this 
was not pursued as large improvements in 
model accuracy were not anticipated by 
enforcing this constraint. Ai'-^.r' mally. 
initial identification tests .'i ed the 
eigenvalues of the L' mat>"ix be neyond 
the frequency range of interest. There- 
fore, the D matrix was cor«Lraineci to be 
D = 1. 


The noise covariance, Q, is also given in 
Fig. 1. Since 0 did not change sub- 
stantially from model 1 to model ?, sub- 
sequent models were identified using the 
0 of model 2. Periodic updates of the 
value of 0 in subsequent models showed no 
large variation in Q. 

Model 1 of Table 1 was used to predict 
engine behavior based upon actual closed 
loop engine test data. This is described 
in the next section. 


~^cs t data Appl ication 

The FlOO engine was tested in the Lewis 
Research Center altitude test facilty to 
evaluate the FlOO Multivariable Control 
(MVC) law (Refs. ? and 3) . During the same 
test period the "Bill of Material" (BOM) 
control was also evaluated as a baseline/ 
backup control model. Thus, there are a 
variety of closed loop operating records 
obtained throughout the flight envelope 
with a number of different power input 
requests. The two multiva-iable data sets 
used in this report were recorded at an 
ALT • 10 000 ft. MN « 0.9 condition as the 
power request was varied (step change) in a 
small (hopefully linear) range about inter- 
mediate engine power. One set corresponds 
to an MVC control test, the other to a BOM 
test. Data were sampled at T . 0.05 sec 
for a 10-second transient, which yields 
K - 200 points for each record in the data 
sets. 



fixed at a given operating point, strong 
system ident if iabi Ity can be guaranteed if 
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Or in other words, if 


det[L] ^ 0 (18) 

The BOM and MVC reference point schedules 
do exhibit the characteristic of (18) , 
therefore a direct identification approach, 
such as IV/AML, can be successfully applied 
to the closed loop input/output data sets 
recorded in the Lewis test facility. 

Sensor 'nstrumentet ion for the input and 
output variables of interest is summarized 
in Table 2. m each case sensor dynamics 
are beyond the 20 radians/sec frequency 
range. Thus, sensor dynamics were ini- 
tially ignored in the identification 
tests. Ambient rcise statistics were ob- 
tained during .Teady-statt engine opera- 
tion. Standard deviatiofl^ were calculated 
at the operating point for the sensed 
values. Signal to noise ratios (SNR's) 
were estimated based upon these ambient 
noise levels i>nd the deviations of the 
various signals from their operating point 
values. These SNR's are included in Table 
2. In each case the SNR's show the level 
of noise tc be small relative to the 
signal. Therefore, the identification 
results should be guite consistent and 
accurate. 

Normalized WF from the BOM and MVC control 
tests is shown in Fig. 3. This is typical 
of the engine inputs in these tests. Power 
spectrum analysis of these inputs shows a 
slightly higher freguency component in the 
MVC inputs, although more total power is 
contained in the BOM inputs. However, for 
both the BOM and MVC inputs most of the 
power is concentrated below 6 raoia’ s/sec. 

The control inputs of Fig. 3 were used in 
conjunction with the identified model 1 to 
predict engine output. Comparing the pre- 
dicted outputs of model 1 with the actual 
outputs, it was apparent that mooel 1 was 
unacceptable. No output was predicted well 
for either BOM or MVC data. Figure 4 is 
typical of the comparison. Slight dis- 
crepancies between simulation and test data 
cannot account for large discrepancies 
between predicted and actual outputs. 


The BOM and MVC control structures, line- 
arized at an operating point, correspond to 
the structure of Fig. 2. The reference 
point and control blocks are different how- 
ever for the two controls. The structure 
of Fig. 2 Is exactly the structure given in 
Ref. 1. Since each control structure Is 


To investigate this further the IV/AML 
method was applied directly to the closed 
loop test data producing models 2 and 3. 
Model 1 was used as a starting point. As 
illustrated in Fig. 5. model 3 accurately 
reproduces the data, from which it was 
generated (BOM). Model 2 results are simi- 
lar. In fact. Fig. 1 shows the error of 
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a11 the outputs for models 1, 2 and 3 to be 
less than 1 percent. However, comparing 
parameters for models 1, 2 and 3 (see 
Table 1) 1t can be seen that while Aj 
remains essentially unchanged elements of 
Li do change substantially. This Im- 
plies a slightly overparameterized model 
structure which does account for the In- 
ability of model 1 to predict BOM and MVC 
engine data. To determine which elements 
should b" “imlnated to remove the over- 
parameter izat I "r., the following procedure 
was adopted. 

First, a reasonably accurate Initial model 
Is assumed. In this case, models 2 and 3 
were used. An Initial covariance matrix 

Pp - tliag(Pj,^) (19) 

is chosen where the pp^ are small to 
Indicate small uncertainty In the model 
parameters. In the engine example 

Pp - 10'^ (20) 

was used. Next, the IV/AM method is ap- 
plied for only a single iteration to data 
for which the model Is overparameterized. 

The method now will be most sensitive to 
removing the uncertainty Inherent In the 
extra parameters. Now, the diagonal ele- 
ments of Pi which correspond to accu- 
rate parameters will not change. However, 
the diagonal elements of Pi, pn, which 
correspond to extra parameters will change 
significantly. Thus, If 

Pll - Poll " ‘Pol 

where t Is a positive threshold, the 
corresponding parameter ai can be set 
to zero. The threshold was selected as 

c - 0.05 (22) 

for the engine data. 

Three elements of Pi satisfied (21) for 
both the MVC ana BOM data. The correspond- 
ing parameters were eliminated and this new 
structure applied to the simulation data 
generated by the Inputs of (14) and (15). 

The resultant IV/AML Identified model Is 
given as model 4 In Table 1. A comparison 
of average fit error Is quite comparable to 
model 1 with full Bi and, in fact, 
shows Improvement In the PT6 comparison. 

Note that the eigenvalue associated with 
PT6 in model 4 represents a frequency of 
approximately 25 radians/sec which Is 
slightly greater tnan the 20 radians/sec 
natural frequency of the PT6 sensor. 
Obviously, the PT6 sensor dynamics can no 
longer be competely ignored In the Inter- 
pretation of the results. Additionally, 
since one mode models the sensor dynamics, 
a second mode may be required to model the 


PT6 engine mode. This was not pursued at 
this time however. 

When used to predict BOM and MVC output 
data, model 4 was still unsatisfactory. 
Model 4 did predict Nl(MVC), N2(MVC), and 
N2(B0M). However, Nl(BOM) and especially 
PT6 for both data sets were not predicted 
well. The error In PT6 Is somewhat ex- 
pected from sensor and Input bandwidth 
considerations. The Nl(BOM) error was not 
expected however. Figure 6 compares pre- 
dicted N1 data using model 4 to actual 
closed loop Nl(BOM) data. Model 4 pre- 
dicted N1 grossly follows the trend of the 
simulated data. Thus, It appears that the 
dynamic portion of model 4 Is correct. 
However, there must then be large dis- 
crepancies In some of the model 4 gain 
terms. These discrepancies are somewhat 
perplexing since model 4 predicted Nl(MVC) 
but not Nl(BOM). 

Recall, however, that the BOM Inputs are 
larger in magnitude than the MVC Inputs, 
and that model 4 represents linearized 
dynamics. Thus, some nonlinear effects may 
be Inherent in the BOM data. This explana- 
tion is not entirely satisfactory since 
N2(B0M) and N2(MVC) were both predicted. 
Further work to resolve this problem Is 
required. The IV/AML Identification method 
was again utilized to further refine the 
model parameters for the structure of model 
4 using the two sets of experimental closed 
loop data. The purpose of this final 
Iteration Is to Identify a single model 
that can accurately predict both sets of 
engine test data and, hopefully, simulation 
data as well. 

Again model 4 was used as an Initial condi- 
tion In the IV/AML method applied to the 
BOM and MVC data. Models 5 and 6 of 
Table 1 resulted. Both models 5 and 6 fit 
their respective data sets quite well. 
Figures 6 to 8. for example, show a good 
fit of the BOM data by outputs predicted 
using model 5. Similar comparisons to MVC 
data were obtained using model 6. More 
Importantly, when the BOM model 5 Is used 
to predict the MVC data, the comparison 
given In Figs. 9 to 11 Is quite reason- 
able. Thus, model 5 (or equlvilently 
model 6) represents a model which predicts 
a class of Inputs and can be used with con- 
fidence In a control design procedure. 

CONCLUSIONS 

The IV/AML method was applied to both open 
loop simulation and closed loop test data 
of an FlOO turbofan engine. The method 
accurately and consistently Identified 
models from both the simulation and test 
data. Due to the structure of the BOM and 
MVC control laws, the engine model Is 
strongly system Identifiable and con- 
sequently a direct Identification approach 
was used on the closed loop data. 
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A third order model structure was derived 
and found to be overparameterlzeo. Three 
parameters were eliminated by sensitivity 
considerations. The simplified structure 
was found acceptable for fitting both simu- 
lation and test data. Test model accuracy 
Is limited to 6 radians/sec since spectral 
analysis of the inputs shows limited signal 
strength above this frequency. 

Also, Identification results Indicate that 
PT6 related model parameter accuracy Is 
further limited by sensor bandwidth and 
that an additional dynamic mode Is required 
to faithfully model PT6. 

Comparisons showed that models identified 
from simulated data generally predicted 
Nl(MVC), N2(MVC). and N2(B0M) test response 
aoequately. However, predictions of 
PT6(MVC) and PT6(B0M were poor and Nl(BOM) 
showed some discrepancies in dynamics. The 
PT6 differences are attributed to the low 
frequency content of the test Input signals 
(<6 radians/sec) , the bandwidtn of the 
sensor, and the high frequency nature of 
the PT6 mode. However, the difference In 
Nl Is attributable to a difference in simu- 
lated versus actual engine performance. 

This conclusion is accurately portrayed In 
a comparison of Identified models. 

Finally, a simplified model determined from 
BOM data accurately predicted not only BOM 
but also MVC test response data. This 
ability to predict engine performance for 
a class of Inputs generates confidence In 
controls designated from this model. Thus, 

It Is concluded that useful dynamic engine 
models can be obtained from closed 1c ip 
test data using the IV/AML identification 
method. This Identification technique, 
then, represents the first step In an 
automated engine control design process. 
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Figure L - Itentive IV /AML version strudurei 
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Figure i - Compitlson ot actual BOM output with a sim- 
ulated BOM output and an output predicted by modef 1 
obtained from simulation data - unconstrained B^. 
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Figure J. * Comptrlson o( actual BOM output with out- 
put predicted ^ model 3 obtained from BOM data - 
unconstrained Bj. 
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Figure 6. * Comparison of actual BOM output with a 
simulated BOM output and outputs predicted by 
models land 3. 


9 Id 



Figure 7. - Comparison til actual BOM output with out* 
put pradictad ^ modal 5 obtalnad from BOM data ■ 
constrainad Bj. 



Figure 0 . ■ Comparison of actual BOM output with out 
put pradictad by modal $ obtainad from BOM data - 
constrainad B|. 



Figure 9. ■ Comparison of actual MVC output with out- 
put predicted ^ model 5 obtained from BOM dita - 
constrained Bj. 



Figure 1C. - Comparison of actual MVC output with out- 
put predicted by model $ obtained from BOM data - 
constrained B|. 
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Figure 11. - Comperison of actual MVC output with out- 
put predicted by model 5 obtained from BOM data - 
constrained Bj. 


